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Executive  Summary 


The  principal  goal  of  this  research  is  to  advance  the  technology  of  structural  health  management 
for  composite  structures.  The  use  of  composites  for  aerospace  structures  is  increasing  rapidly; 
however,  composite  structures  are  susceptible  to  impact  damage  and  delaminations  and  cracks 
may  reach  critical  length  before  visual  detection.  Wave  propagation  based  methods  have  shown 
promise  for  SHM  of  composite  structures.  The  spectral  finite  elements  (SFE)  method  is  highly 
suitable  for  wave  propagation  analysis  due  to  its  frequency  domain  approach  which  yields 
models  that  are  many  orders  smaller  than  conventional  FEM.  Also,  the  frequency  domain 
formulation  of  WSFE  enables  direct  relationship  between  output  and  input  through  system 
transfer  function.  Wavelet  spectral  finite  element  (WSFE)  method  overcomes  the  “wrap  around” 
problem  to  accurately  model  2-D  plate  structures  of  finite  dimensions,  unlike  the  existing  Fourier 
transform  based  SFE.  The  specific  accomplishments  of  the  research  are: 

•  Developed  WSFE  model  for  damaged  composite  plate  with  transverse  crack; 

•  Validated  WSFE  modeling  of  Lamb  wave  propagation  in  healthy  composite  plates  through 
experimental  measurements; 

•  Implemented  ‘baseline-free’  Damage  Force  Indicator  method  for  delamination  detection 
using  dynamic  stiffness  matrix  from  WSFE  model  of  healthy  plate; 

•  Developed  Modified  Time  Reversal  method  for  Lamb  wave  based  ‘baseline-free’  damage 
diagnostics;  and 

•  Investigated  effect  of  tone  burst  center  frequency  on  instantaneous  phase  based  delamination 
detection. 

This  award  partially  supported  2  MS  students  and  1  undergraduate  honors  student.  The  research 
resulted  in  the  publication  of  6  conference  papers  and  4  journal  articles  (in  review/preparation). 


1.  BACKGROUND  AND  SIGNIFICANCE 


The  use  of  composites  for  aerospace  structures  is  increasing  rapidly  due  to  several  advantages 
such  as  lighter  weight,  fewer  joints,  improved  fatigue  life,  and  higher  resistance  to  corrosion. 
However,  composite  structures  have  several  forms  of  damage  such  as  delamination  between 
plies,  fiber-matrix  debonding,  fiber  breakage,  and  matrix  cracking.  These  damages  often  occur 
below  the  surface  due  to  fatigue,  foreign  object  impact,  etc.,  and  may  not  be  visible.  The  existing 
paradigm  for  structural  safety  leads  to  expensive  inspections,  unnecessary  downtime  and 
retirement,  and  sometimes  catastrophic  failures  without  any  warning.  Therefore,  a  structural 
health  management  (SHM)  system  capable  of  performing  both  diagnostics  and  prognostics  for 
composite  structures  is  urgently  needed.  A  reliable  SHM  system  will  provide  tremendous 
benefits  in  terms  of  life-cycle  costs  and  allow  realization  of  CBM+. 

Lamb  waves  (fundamental  symmetric  and  asymmetric  modes)  have  shown  promise  for  damage 
diagnosis  of  composite  structures.  The  fundamental  idea  behind  Lamb  wave  propagation  based 
diagnostics  is  that  different  types  of  damages  interact  differently  with  waves.  Therefore,  based 
on  the  measured  time  history  of  the  propagated  wave,  the  traveling  time,  speed  reduction,  and 
wave  attenuation  parameters  are  extracted  and  used  as  the  damage  identification  variables.  Due 
to  the  dispersive  nature  of  the  Lamb  waves,  accurate  identification  of  localized  events  is  required 
to  determine  the  location  of  the  damage  in  the  structure.  Advanced  signal  processing  using 
wavelet  transform,  Hilbert-Huang  transform,  etc.  helps  damage  recognition. 

Conventional  finite  element  method  (FEM)  is  not  suited  for  wave  propagation  problems  since 
the  element  size  should  be  comparable  to  wavelength  which  is  very  small  at  high  frequencies. 
This  results  in  large  system  size  and  enormous  computational  cost.  In  addition,  solving  inverse 
problems  (as  required  for  diagnostics)  is  very  difficult  using  FEM.  Spectral  finite  element  (SFE), 
which  follows  FE  modeling  procedure  in  the  transformed  frequency  domain,  is  highly  suitable 
for  wave  propagation  analysis.  SFE  models  are  many  orders  smaller  than  FEM  and  especially 
suitable  for  on-board  diagnostics.  The  frequency  domain  formulation  of  SFE  enables  direct 
relationship  between  output  and  input  through  system  transfer  function  (or  FRF). 

Wave  propagation  analysis  in  composite  plates  using  wavelet  spectral  finite  element  has  been 
investigated  in  this  research.  The  2-D  Wavelet  based  SFE  (WSFE)  can  model  2-D  plate 
structures  of  finite  dimensions  accurately.  This  is  due  to  the  use  of  localized  Daubechies  scaling 
functions  as  basis  for  approximation  of  the  spatial  dimension.  Presence  of  elastic  coupling  in 
anisotropic  materials  results  in  coupled  governing  differential  equations.  The  final  reduced  ODEs 
after  spatial  and  temporal  approximations  are  in  the  form  of  a  set  of  coupled  ODEs,  for  each 
temporal  and  spatial  sampling  point.  The  solution  of  these  ODEs  to  derive  the  exact  shape 
function  involves  determination  of  wave  numbers  and  the  amplitude  ratio  matrix.  Unlike 
isotropic  cases,  the  process  of  solution  for  composite  materials  is  more  complicated  and  is  done 
by  posing  it  as  polynomial  eigenvalue  problem  (PEP). 
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2.  Wavelet  Spectral  Finite  Element  Modeling  with  Transverse  Crack 


The  governing  differential  equations  of  a  homogeneous,  symmetric  and  balanced  laminate  based 
on  first  order  shear  deformation  theory  (FSDT)  can  be  expressed  in  terms  of  displacements 
(u  ,v  ,w  ,0 ,(//  in  Fig.  1)  as, 
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Fig.  1 .  The  plate  element  with  associated  DOF 
The  associated  boundary  conditions  are, 


Nxx=Nnx+N„nv, 


Nyy=N3Vnx+Nxyny, 


Vx=Vxnx+Vyny 


Mxx=-MxxnX~Mxyny’ 


Myy=-Mxynx~Myyny 


(6) 

(7) 
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Where  Nxx  and  N  are  the  applied  normal  forces  in  X  and  Y  direction,  M  xx  and  M  are  the 
applied  moments  about  Y  and  X  axes  and  Vx  is  the  applied  shear  force  in  Z  direction.  The 
terms  nx  and  n  are  unit  normal  into  X  and  Y  directions,  respectively. 


An  assumption  is  made  that  the  lamina  is  specially  orthotropic  which  is  true  for  many 
commercially  available  unidirectional  laminae.  Transverse  motion  is  studied  and  therefore  the  in¬ 
plane  displacements  can  be  left  out  in  the  formulation  hence  forth.  Further  discussion  of  the 
WSFE  formulation  can  be  found  in  [1,  4-6]  and  is  not  reproduced  here.  Only  the  final 
transformed  governing  equations  and  boundary  conditions  are  mentioned  for  completeness.  The 
transformed  governing  equations  that  have  transverse  displacement  components  are, 
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and  the  boundary  conditions  can  be  written  as, 
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(13) 


where  m  and  n  number  of  discretization  of  spatial  and  time  windows. 

The  schematic  diagram  of  spectral  plate  element  (SPE)  with  a  transverse  open  and  non¬ 
propagating  crack  is  shown  in  Fig.  2.  The  length  of  the  element  in  the  X  direction  is  Lx  and  the 

plate  has  finite  dimension  of  Ly  in  the  Y  direction.  The  crack  is  located  at  a  distance  of  Lx  from 
the  left  edge  of  the  plate  and  has  a  length  2c  in  the  Y  direction. 
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Fig.  2.  The  plate  element  with  transverse  crack 

Two  sets  of  nodal  displacement  fields,  { wl ,  <f\ ,  \j/x  }and  {w2,(j)2,y/2},  are  assumed  on  the  left  and 
right  hand  side  of  the  crack  whose  wavenumber  domain  representations  are  { w, ,  ^  ,\frx  }and 
{w2,(j)2,y/2}  and  they  can  be  expressed  as, 
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matrix  with  diagonal  terms 
[f?]is  3X6  amplitude  ratio  matrix  and  a  total  of 
12  unknown  constants  from  { C, }  and{C2}  to  be  determined  using  boundary  conditions. 

The  unknown  constants  [  C,  j  and  {C2}  can  be  calculated  as  a  function  of  the  nodal  spectral 
displacements  using  the  boundary  conditions  as  follows. 

1 .  At  the  left  edge  of  the  element  (x  =  0) 

w,  =<?p  A=q2,  (15) 

2.  At  the  crack  location  (  x  =  Lx  for  w,  and  x  =  0  for  w2 ) 

w,  =w2,  i//,  =  y/2  (Continuity  of  transverse  displacement  and  slope  across  crack) 

$  -  <fi2  =  OMxx  (Discontinuity  of  slope  across  crack) 

Mxxl  =  M xx2,  M  j  =  Mxy 2  (Continuity  of  bending  moments) 

Vxl  =  Vx2  (Continuity  of  shear  force)  (16) 
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3 .  at  the  right  edge  of  the  element  ( x  =  L  -  Lx  for  w2 ) 

w2=q4’  &=q5>  ¥i=q6  (17) 

6  is  the  flexibility  of  the  plate  along  the  crack  edge  and  the  way  of  calculating  it  is  presented  in 
the  next  section.  MxxA  ,  Mxrt  and  Vxl  are  the  moments  and  shear  force  obtained  in  terms  of 

w2 ,  (j)2 ,  i//2  and  similar  for  M xx2 ,  M  2  and  Vx2 .  These  boundary  conditions  can  be  further 
expressed  in  matrix  form  as, 
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where  ux  =  {qt,q2,q3}T  ,  u2  =  [q4,q5,q6}T  ,M3,M4  eC3X6and  M2,M3  eC“6 

The  shear  force  and  bending  moments  at  the  left  and  right  boundaries  of  the  plate  can  be  written 
as, 
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Combining  eqns.  (18)  and  (19), 

=[K} 


u , 


Ur , 


(20) 


K  is  the  frequency  dependent  dynamic  stiffness  matrix  of  the  SPE  with  transverse  open  and  non¬ 
propagating  crack. 

The  bending  flexibility  on  both  sides  of  the  crack  can  be  expressed  in  dimensionless  form. 

albf  (t)^(  t),  y=j-  (21) 

Ly 

where  His  the  total  thickness  of  the  plate,  F(y)  is  a  correction  function  and  a°bb  is  the 
dimensionless  bending  compliance  coefficient  at  the  crack  center  which  can  be  expressed  as, 
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for  the  range,  0  <  £  <  0.7 .  Here  h  is  the  depth  of  crack  at  any  given  y  and  \  is  the  central  crack 
depth.  / (y)  represents  the  crack  shape  function  defined  as, 
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where  e  is  the  base  of  natural  logarithm.  The  correction  factor  is  given  as, 

=  2c///  +  3(v  +  3)(l-v)4[l-/(y)] 

Vy  2c/JH'  +  3(v  +  3)(l-v)a^ 


Fig.  (4)  presents  the  variation  of  6  for  few  different  values  of  crack  parameters. 


(23) 


(24) 


Fig.  4.  Variation  of  6  withy  along  the  crack 


The  formulated  2-Element  WSFE  model  with  a  transverse  crack  is  used  to  study  the  wave 
scattering  in  graphite-epoxy  AS4/3501  composite  plate  (Fig.  5).  Time  domain  results  (Fig.  6) 
show  partial  reflections  of  the  excited  wave  due  to  the  crack  at  different  time  with  varying 
distances  to  the  crack  from  the  excitation  point.  The  presence  of  crack  does  not  change  the  peak 
amplitude  and  the  group  velocity  at  which  the  pulse  is  propagated.  As  expected,  an  additional 
reflection  from  the  crack  arrives  before  the  reflection  from  the  fixed  boundary.  The  wave 
propagation  in  space  at  different  time  instants  (with  crack  at  Lx  =  0.50  m  away  from  the 

excitation  point)  shows  partial  wave  reflection  from  crack  (Figure  7).  Experimental  validation  of 
these  results  is  currently  under  progress. 
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Fig.  5  Composite  plate  with  transverse  crack 


Time  (m  s) 

Fig.  6  Transverse  velocity  for  various  crack  locations 


Fig.  7  Wave  propagation  (tone  burst)  in  composite  plate  with  crack  at  x=0.5  m  at  630ps  (left)  and  750  ps 
(right) 


3.  Composite  Delamination  Detection  Using  WSFE  and  Damage  Force  Indicator  Method 

The  Damage  Force  Indicator  (DFI)  Method  was  implemented  for  detection  of  delamination 
without  the  need  for  baseline  data.  The  basic  idea  in  the  DFI  method  is  that,  using  the  dynamic 
stiffness  of  healthy  structure  and  the  measured  damage  response,  the  force  necessary  to  create 
damage,  called  the  ‘damage  force,’  can  be  computed.  These  damage  forces  are  then  normalized 
and  plotted  at  all  sensor  locations.  The  damage  force  indicator  peaks  only  at  those  locations 
which  are  adjacent  to  the  damage,  which  clearly  helps  to  locate  the  damage.  Some  of  the  major 
advantage  of  the  DFI  method  are  that  it  does  not  need  any  baseline  measurement,  knowledge  of 
input  is  not  necessary  and  it  is  computationally  very  fast,  making  it  an  ideal  tool  for  on-line 
health  monitoring  when  the  dynamic  stiffness  matrix  is  obtained  from  WSFE  model. 
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The  experimental  setup  (Fig.  8)  and  results  for  Lamb  wave  measurements  are  discussed  in  this 
section.  The  composite  plate  used  in  this  study  is  fabricated  from  AS4/3501-6  carbon-epoxy  pre- 
preg  using  vacuum  bagging  and  oven  curing  technique.  The  plate  has  dimensions  of  254  x  127  x 
1.25  mm  (10  x  5  x  0.05  in)  and  consists  of  8  cross-ply  [0/90]2S-  A  piezoelectric  (PZT)  actuator 
(diameter  13.5  mm  and  thickness  0.22  mm)  is  affixed  onto  the  composite  plate  using  epoxy.  A 
National  Instruments  PXI  6339  and  a  BNC-2110  board  are  used  to  generate  signals  and  a 
QuickPack®  power  amplifier  is  used  to  amplify  the  actuation  signal.  Tone  burst  actuates  at  15 
kHz  with  3.5  cycles  to  generate  Ao  Lamb  waves  (Fig.  9).  A  Scanning  Laser  Doppler  Vibrometer 
(SLV)  is  employed  to  acquire  the  signals  at  designated  locations.  Based  on  Lamb  wave 
dispersion  curve,  phase  velocity  calculations  indicate  that  the  acquired  signals  are  Lamb  wave  Ao 
mode. 


Fig.  8.  Experimental  setup  (L)  and  composite  plate  with  attached  PZT  actuator  (R) 


Frequency  (kHz) 


Fig.  9.  15  kHz  3.5  cycle  tone  burst  in  time  domain  (left)  and  its  frequency  domain  representation 

(right) 
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Experimental  data  was  compared  with  WSFE  predictions  (Fig.  10)  for  a  healthy  composite  plate 
to  validate  WSFE  modeling.  Damage  is  generally  inferred  by  comparing  healthy  and  damaged 
(delaminated)  responses  from  a  structure.  At  a  sensor  location  close  to  delamination,  there  is  a 
clear  difference  in  the  second  wave  packet  (Fig.  11).  For  DFI  method,  we  use  only  dynamic 
stiffness  matrix  from  WSFE  model  (not  healthy  plate  response  data)  and  delaminated  plate 
response.  The  DFI  method  correctly  identifies  damage  since  the  dominant  peak  is  located  at 
degree  of  freedom  15  which  belongs  to  the  sensor  located  right  after  the  delamination  (Fig.  12). 
Further  details  are  given  in  [1,  7-8], 


Fig.  10.  Comparison  of  velocity  response  from  WSFE  and  experiments  at  two  sensor  locations 


Fig.  11.  Velocity  responses  from  healthy  and  delaminated  plates  (sensor  location  close  to 

delamination) 
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Fig.  12.  Damage  force  indicator  for  composite  plate  (a)  for  healthy  (b)  damaged  plate 


4.  Modified  Time  Reversal  method  for  Lamb  wave  based  ‘baseline-free’  damage 
diagnostics 

Time  reversal  method  (TRM)  is  a  new  approach  to  mitigate  Lamb  wave  dispersion  effects  and 
increase  their  applicability  for  SHM.  We  have  developed  modified  time  reversal  method  (Fig. 
13),  wherein  some  of  the  steps  are  modified  (as  denoted  by  bold  font)  compared  to  TRM. 

1 .  A  signal  is  sent  from  transducer  A  and  recorded  at  transducer  B 

2.  The  received  signal  at  B  is  reversed  in  time,  that  is  vb$)  -*■  Vb(-«) 

3.  The  time  reversed  signal  is  sent  from  A  to  B  where  it  is  recorded 

4.  The  received  signal  at  B  is  time  reversed  and  compared  to  the  original  signal  actuated 


6)  Compare  signals 


Figure  13:  Modified  time  reversal  method  using  pitch-catch  arrangement  of  transducers 

Based  on  the  assumption  of  linear  reciprocity  in  the  time  reversal  method,  the  signal  due  to  the 
time  reversal  process  and  the  original  signal  compare  exactly;  however,  for  nonlinearities 
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introduced  into  the  system  (damage),  the  linear  reciprocity  of  the  system  breaks  down  and  the 
irregularities  between  the  two  signals  indicate  the  presence  of  damage.  Unlike  most  other  health 
monitoring  techniques,  this  time  reversal  process  introduces  a  baseline-free  method  for  damage 
detection.  In  our  MTRM,  a  major  benefit  is  that  only  one  actuator  is  necessary  instead  of  two  for 
any  actuator/sensor  pair.  We  have  shown  mathematically  that  the  final  result  is  same  as  the  TRM 
process  adopted  by  other  researchers.  Figure  14  shows  the  experimental  set-up  for  composite 
plate  with  impact  damages  (left)  and  MTRM  results  for  various  damage  levels.  Results  show  that 
magnitudes  of  both  damage  indices  are  directly  correlated  to  the  severity  of  damage.  Further 
details  are  given  in  [3,  9-10]. 
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Figure  14.  Experimental  set-up  for  composite  plate  with  impact  damages  (left)  and  effect  of 
impact  on  damage  index  (right).  Consistent  increase  in  both  damage  indices  with  increasing 
damage  is  observed. 

5.  Effect  of  tone  burst  center  frequency  on  instantaneous  phase  based  delamination 
detection 

Lamb  waves  in  plates,  known  as  guided  waves,  are  elastic  perturbations  that  can  propagate  for 
long  distances  in  thin  plate-like  structures  with  little  amplitude  loss.  But  Lamb  waves  are  highly 
dispersive,  in  other  words  wave  propagation  depends  on  the  frequency-thickness  relationship. 
Therefore,  tone  burst  signals  wherein  most  of  the  energy  is  concentrated  at  a  central  frequency, 
are  generally  used  for  Lamb  wave  excitation.  Several  parameters  such  as  central  frequency, 
number  of  cycles  and  amplitude  of  the  incident  Lamb  wave  may  have  significant  impact  on 
damage  detection.  Published  papers  discuss  the  effects  of  actuation  pulse  parameters  in  a  general 
way,  but  specifics  are  seldom  reported.  In  our  research,  frequency  of  tone  burst  signals  is  varied 
from  20  to  80  kHz  with  5  kHz  increments  with  fixed  number  of  cycles  (3.5  cycles)  and  actuator 
voltage  input.  Piezoelectric  (PZT)  transducers  are  used  as  actuators  and  both  PZT  and  LDV 
perform  sensing.  LDV  collects  responses  of  flexural  waves  and  PZTs  collect  responses  of  both 
flexural  waves  and  axial  waves.  Those  responses  are  evaluated  with  three  different  methods, 
namely  comparisons  of  energy,  frequency  and  phase  between  healthy  and  delaminated  plates.  A 
10  mm  wide  delamination  is  created  at  the  mid-plane  along  the  width  of  8-layer  cross-ply  plate. 
Wavelet  transform  (WT)  is  used  for  frequency  evaluation  and  both  WT  and  Hilbert-Huang 
transform  (HHT)  are  used  for  phase  comparisons. 
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The  composite  plate  used  in  this  study  has  dimensions  of  432  x  324  x  1.25  mm  which  is 
fabricated  from  8  plies  of  AS4/3501-6  pre-preg  using  vacuum  bagging  and  oven  curing 
technique.  Hi-temp  mold  release  wax  from  PARTALL  was  applied  between  4th  and  5th  layers  to 
create  the  delamination  area  (10  x  152  mm).  Six  circular  PZTs  (7  x  0.2  mm)  are  affixed  onto  the 
composite  plate  at  50  mm  distance  from  each  other  using  epoxy  (Fig.  15).  An  Agilent  33220A  20 
MHz  function/ Arbitrary  waveform  generator  is  used  to  generate  tone  burst  signals  and  an  AA 
Lab  systems  piezo-amplifier  is  used  to  amplify  the  excitation  signals.  Tone  burst  excitation 
signals  (3.5  cycles)  are  sent  to  each  PZT  at  13  different  frequencies  from  20  to  80  kHz  in  5  kHz 
increments.  PZT  sensors  are  connected  to  BK  Precision  2540  60  MHz  digital  oscilloscope  which 
sends  stored  signals  to  computer. 


432  mm 
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Figure  15.  AS4/3501-6  pre-preg  cross  ply  composite  plate  with  PZTs  (left)  and  schematic 

diagram  with  delamination  area  (right) 

In  HHT,  empirical  mode  decomposition  (EMD)  is  used  to  decompose  original  signals  into  N 
intrinsic  mode  functions  (IMFs).  The  Hilbert  transform  of  individual  IMF’s  gives  instantaneous 
phase  and  the  total  instantaneous  phase  is  obtained  by  a  sum  of  the  instantaneous  phases  for  each 
IMFs.  Phase  difference  calculation  with  HHT  shows  that  the  signal  for  Path  2  has  significant 
phase  change  compared  to  other  paths  (Fig.  16).  To  quantify  the  magnitude  of  phase  change, 
Figure  17  shows  the  mean  of  phase  difference  values  for  each  path  (with  Path5  as  reference). 
Path2  (delamination  area)  shows  much  larger  phase  difference  and  indicates  damage  very 
clearly.  The  phase  difference  is  somewhat  smaller  at  lower  frequencies.  Thus  the  experimental 
results  for  the  selected  frequency  range  (20  kHz  to  80  kHz  with  5  kHz  increments)  show 
significant  effects  of  actuation  frequency  on  both  energy  and  phase  difference  angle  [2,  11-12], 
Further  studies  will  include  understanding  the  relationship  between  actuation  frequencies  and 
energy  and  phase  difference  observed  and  different  size  of  delamination  area. 
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Fig.  16.  Phase  difference  calculation  with  HHT 


■  Path  1  vs.  Path  5 

■  Path  2  vs.  Path  5 

■  Path  3  vs.  Path  5 

■  Path  4  vs.  Path  5 
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Figure  17.  Phase  angle  difference  versus  frequency 
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